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The mechanism of non- Abelian color confinement is studied in SU(2) lattice gauge theory in terms 
of the Abelian fields and monopoles extracted from non- Abelian link variables without adopting gauge 
fixing. Firstly, the static quark-antiquark potential and force are computed with the Abelian and 
monopole Polyakov loop correlators, and the resulting string tensions are found to be identical to 
the non-Abelian string tension. These potentials also show the scaling behavior with respect to 
the change of lattice spacing. Secondly, the profile of the color-electric field between a quark and 
an antiquark is investigated with the Abelian and monopole Wilson loops. The color-electric field 
is squeezed into a flux tube due to monopole supercurrent with the same Abelian color direction. 
The parameters corresponding to the penetration and coherence lengths show the scaling behavior, 
and the ratio of these lengths, i.e, the Ginzburg-Landau parameter, indicates that the vacuum type 
is near the border of the type 1 and type 2 (dual) superconductor. These results are summarized 
that the Abelian fundamental charge defined in an arbitrary color direction is confined inside a 
hadronic state by the dual Meissner effect. As the color-neutral state in any Abelian color direction 
corresponds to the physical color-singlet state, this effect explains non-Abelian color confinement 
and supports the existence of a gauge-invariant mechanism of color confinement due to the dual 
Meissner effect caused by Abelian monopoles. 

PACS numbers: 12.38.AW,14.80.Hv 



I. INTRODUCTION 

Color confinement in quantum chromodynamics 
(QCD) is still an important unsolved problem [lj. 
't Hooft and Mandelstam [1] conjectured that the 
QCD vacuum is a kind of a dual superconducting state 
caused by condensation of magnetic monopoles. The 
color charges are then confined inside hadrons due to 
formation of the color-electric flux tube through the dual 
Meissner effect. However, in contrast to the Georgi- 
Glashow model or SUSY QCD [§[ with scalar fields, 
it is not straightforward to identify the color-magnetic 
monopoles in QCD. 

An interesting idea to realize this conjecture is pro- 
posed by 't Hooft 0], such that SU(3) QCD can be re- 
duced to an Abelian [U(l)] 2 theory by adopting a par- 
tial gauge fixing, and the color-magnetic monopoles ap- 
pear according to tt 2 (SU(3)/[U(1)] 2 ) = Z 2 . The role of 
monopoles for the confinement mechanism is investigated 
extensively on the lattice by applying Abelian projec- 
tion in the maximally Abelian (MA) gauge [1, 0| , where 
monopoles are extracted a la DeGrand-Toussaint fl~oj ] as 
in compact U(l) lattice gauge theory. It is then found 
that the results strongl y su pport the dual superconduct- 
ing scenario [111, d El, Hi, EE EE E3, El ■ The confining 
properties are dominated by the Abelian fields [HI, EE Elf 
and monopoles [TJ, HE SI S3, El, H , which are called 
Abelian dominance and monopole dominance, respec- 
tively. The color-electric flux is squeezed by the dual 
Meissner effect [H, EE E3, El, EH • Moreover monopole 



condensation is confirmed by the energy-entropy balance 
of the monopole trajectories [IE HE E3, HI]- These re- 
sults indicate that there must exist a dual Ginzburg- 
Landau t ype theory as an infrared effective theory of 

QCD [H,|3E[3l|. 

However, there are still serious problems to prove this 
scenario. Firstly, there are infinite ways of the partial 
gauge fixing. Since the behavior of the monopoles can 
depend on the gauge choice, it is not clear if the lat- 
tice results in the MA gauge are universal. Note that 
in the Polyakov (PL) gauge, 't Hooft's color-magnetic 
monopoles Q propagate only in the time direction, which 
cannot confine static color charges [H]. Secondly, as 
the 't Hooft scheme essentially uses the Abelian degrees 
of freedom, it is not explained how non-Abelian color 
charges are confined. 

Recently, we have obtained clear numerical evidences 
of Abelian dominance and the dual Meissner effect in 
local unitary gauges such as the F12 and the PL gauges 
in SU(2) lattice gauge theory [33| . where we have used the 
DeGrand-Toussaint monopoles [l0| as in the MA gauge. 
These results provide us with the following ideas. 

1. The DeGrand-Toussaint monopoles on the lat- 
tice can be different from the 't Hooft color- 
magnetic monopoles [7|. 

2. There must exist a gauge-invariant mechanism of 
color confinement due to Abelian monopoles [34l . 

m. 

In this paper, we aim to show detailed numerical evi- 
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dence how these ideas are realized. We investigate the 
confining properties in SU(2) lattice gauge theory in 
terms of the gauge- variant Abelian fields and monopoles 
extracted from non- Abelian link variables without adopt- 
ing any spatially local or non-local gauge fixing. We find a 
gauge-invariant Abelian mechanism of color confinement 
due to Abelian monopoles work even in the continuum 
limit of SU(2) QCD, although we consider gauge- variant 
Abelian operators. The results may also apply to SU(3) 
gauge theory, since the essential features are not altered. 

The paper is organized as follows. In Sec. [Til we ex- 
plain how to extract the Abelian fields and the monopoles 
from non-Abelian link variables without gauge fixing. 
In Sec. HH and Sec. IIV| we compute the static quark- 
antiquark potential and the force with the Abelian and 
monopole Polyakov loop correlators, and find that the 
string tensions exhibit Abelian dominance and monopole 
dominance. These potentials also show the scaling be- 
havior with respect to the change of lattice spacing. 
In Sec. [V] we investigate the correlation function be- 
tween the Abelian operators and the Wilson loop. We 
observe that the color-electric field is squeezed into a 
flux tube due to monopole supercurrent with the same 
Abelian color direction. The parameters corresponding 
to the penetration depth and the coherence length show 
the scaling behavior, and the ratio of these lengths, i.e, 
the Ginzburg-Landau (GL) parameter, indicates that the 
vacuum type is near the border of the type 1 and type 2 
(dual) superconductor. In Sec. IVI1 we discuss implica- 
tions of our results, i.e, the Abelian fundamental charge 
defined in an arbitrary color direction is confined by the 
dual Meissner effect. As the color-neutral state in any 
Abelian color direction corresponds to the physical color- 
singlet state, the dual Meissner effect for the Abelian 
fundamental charge can also explain confinement of non- 
Abelian color charges. The final section IVIII is devoted 
to conclusion and remarks. Our preliminary results are 
already published in Ref. [36j |. 

II. ABELIAN PROJECTION AND 
EXTRACTION OF MONOPOLES 

We explain how to extract the Abelian fields and 
the color-magnetic monopoles from the thermalized non- 
Abelian SU(2) link variables, 

U fl (s) = U°(s)+ia-U IM (s) , (1) 

where a = (a- 1 , a 2 , a 3 ) is the Pauli matrix. Abelian link 
variables in one of the color directions, for example, in 
the a 1 direction are defined as 

Un(s) — cos M (s) + ia 1 sin M (s) , (2) 

where 

^( S )=arctan(^) (3) 



correspond to the Abelian fields. Without gauge fixing 
the Abelian fields in any color directions should be equiv- 
alent. 

We then define the Abelian field strength tensors as 

= Q flv (s)+2m tlu {s) , (4) 

where <E [— tt, n] and n^(s) is an integer correspond- 
ing to the number of the Dirac strings piercing the pla - 
quette. The monopole currents are then defined by [10| 

47T 

= -^upadp,n pa (s + i>) G Z, (5) 
where dp is regarded as a forward difference. 

III. ABELIAN DOMINANCE 

We show the result of the Abelian static potential 36] . 
We generate thermalized gauge configurations using the 
SU(2) Wilson action at a coupling constant (3 = 2.5 on 
the lattice Nf x N t — 24 3 x 24, where the lattice spacing 
a{[3 — 2.5) = 0.0836(8) [fm] is fixed by assuming y/cr = 
440 [MeV]. 

By using the multi- level noise reduction method [3^ |. 
we evaluate the Abelian static potential Va from the cor- 
relation function of the Abelian Polyakov loop operator 

N t -i 

P A = exp[z Hs + kl)] , (6) 

separated at a distance R as 

Va(Ri) = -^-HPa(0)PX(R)) • (7) 
aN t 

The q-q distance R is improved to Rj = {A-kG{R))^ 1 in 
order to reduce the lattice artifact due to finite-lattice 
spacing, where G(R) is the Green function of the lattice 
Laplacian in three dimensions [38l . [3^ |. For the multi- 
level method, the number of sub-lattices adopted is 6 
and the sublattice size is 4a. 

The result is plotted in Fig. [T] where the non-Abelian 
static potential computed from the ordinary Polyakov 
loop correlation function is also plotted for compari- 
son. The number of independent gauge configurations 
is iV COI1 f = 10 in both cases, but the number of inter- 
nal updates in the multi-level method is 15000 for the 
non-Abelian case and 160000 for the Abelian case. The 
statistical errors are determined by the jackknife method. 

We fit the potential to the usual functional form 

V & t(R) = o-R-c/R + n, (8) 

where a denotes the string tension, c the Coulombic co- 
efficient, and fi the constant. The result is summarized 
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FIG. 1: The Abelian static potential in comparison with the 
non-Abelian one. The lines denote the best fitting curve to 
the function Vat(-R)- 



TABLE I: Best fitted values of the string tension aa 2 , the 
Coulombic coefficient c and the constant jj,a. Vna denotes 
the non-Abelian static potential. N[ up is the number of inter- 
nal updates in the multi-level method. FR means the fitting 
range. The x 2 for the central value is xV-^df < O.f . 





aa 2 


c 




FR(R/a) 


Mup 


Vna 
Va 


0.0348(7) 
0.0352(16) 


0.243(6) 
0.231(39) 


0.607(4) 
1.357(17) 


3.92 - 9.97 
4.94 - 9.97 


15000 
160000 



TABLE II: Simulation parameters for the measurement of the 
static potential and the force from Pa, P P h and P mon . A^rgt 
is the number of random gauge transformations. 



p 


x N t 


a(/3) [fm] 


N con { 


•Wrgt 


2.20 


24 3 x 4 


0.211(7) 


6000 


1000 


2.35 


24 3 x 6 


0.137(9) 


4000 


2000 


2.35 


36 3 x 6 


0.137(9) 


5000 


1000 


2.43 


24 3 x 8 


0.1029(4) 


7000 


4000 



IV. MONOPOLE DOMINANCE 

A. The monopole Polyakov loop 

Wc investigate the monopole contribution to the static 
potential in order to examine the role of monopolcs 
for confinement. The monopole part of the Polyakov 
loop operator is extracted as follows. Using the lattice 
Coulomb propagator D(s-s'), which satisfies d v d' v D{s — 
s') = —8 ss i with a forward (backward) difference d v (d' v ), 
the temporal component of the Abelian fields 6*4 (s) are 
written as 

9 i { S ) = -Y J D{ S -s')[d l p Q vi {s l )+d i {d l p e u { s '))] . (11) 

s' 

Inserting Eq. ([11]) (and then Eq. |@J)) to the Abelian 
Polyakov loop (J6j) , we obtain 



in Table Q] We find Abelian dominance such that the 
Abelian string tension is the same as the non-Abelian 
one. 

Here we make a comment on the theoretical observa- 
tions of the Abelian dominance of the string tension us- 
ing the character expansion Hl[ . The authors of [4l[ 
say that they have proved exactly the Abelian dominance 
without gauge- fixing based on a relation for any two half- 
integer representations ji > ji 



km W ^ T \ = 0, 
t^oo W h [R, T] 



(9) 



where 



N t -1 

P ph = exp{-i £ + s')^0-4(s')} - 



k=0 



N t -1 



P mon = exp{-27ri D ( s + kl ~ s ')5>,4(s')} • 

(12) 



fc=0 s' 



We call Pph the photon and P m0 n the monopole parts of 
the Abelian Polyakov loop, respectively [2j|. The latter 
is due to the fact that the Dirac strings n V 4,{s) lead to 
the monopole currents in Eq. (Tcj . Note that the sec- 
ond term of Eq. (fTTj) does not contribute to the Abelian 
Polyakov loop in Eq. ^ . 



WAG] 



2j 



-<Xj[U(C)}> 



(10) 



B. Simulation parameters 



and Xj\jd] * s the SU(2) group character in representation 
j . In Ref . [4l[ , the above relation is derived on some con- 
siderations of screening effects and gluelump energy. But 
nobody knows an exact method of analytic calculations 
in the infrared region of QCD, so that their considera- 
tions about the screening effects and gluelump energy are 
not exact theoretically, although very plausible. Our nu- 
merical observations here are hence non-trivial and they 
suggests the above relation © is actually exact. 



We then compute the static potential from the 
monopole Polyakov loop correlation function. However, 
since Eq. (fl2|) contains the non-local Coulomb propaga- 
tor D(s — s') and the Polyakov loop is not written as a 
product of local operators along the time direction, the 
multi-level method cannot be applied. Without such a 
powerful noise reduction method, it is hard to measure 
the Polyakov loop correlation function at zero tempera- 
ture with the present available computer resource. Thus 
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FIG. 2: The static potential (top) and the force (bottom) 
from the non-Abelian, the Abelian, the monopole and the 
photon Polyakov loop correlation function at (3 = 2.35 on the 
24 3 x 6 lattice. 



FIG. 3: The same plot as in Fig.[2]at (3 = 2.35 on the 36 3 x 6 
lattice. 



we consider a finite temperature T / system in the 
confinement phase. We set T = 0.8 T c . In order to ex- 
amine the scaling behavior of the potential, we simulate 
the Wilson action on the 24 3 x (N t — 4, 6, 8) lattices. 
We choose the gauge coupling for each N t so as to keep 
the same temperature. We also investigate the spatial 
volume dependence of the potential for the N t = 6 case. 
Simulation parameters are summarized in Table ITT1 The 
lattice spacing a((3) is determined by using the Sommcr 
scale (tq = 0.5 [fm]) at zero temperature. 

C. Noise reduction by gauge averaging 

Since the signal-to-noise ratio of the correlation func- 
tions of Pa, Pph and P mon are still very small with no 
gauge fixing, we adopt a new noise reduction method [36| . 
For a thermalized gauge configuration, we produce many 
gauge copies applying random gauge transformations. 
Then we compute the operator for each copy, and take 
the average over all copies. It should be noted that 
as long as a gauge-invariant operator is evaluated, such 
copies are identical, but they are not if a gauge-variant 
operator is evaluated as in the present case. The results 



obtained with this method are gauge-averaged, thus, 
gauge- invariant . 

In practice, we prepare a few thousand of gauge copies 
for each independent gauge configuration (see Table HT1) . 
We also apply one-step hypercubic blocking (HYP) [42| 
to the temporal links for further noise reduction. The 
short-distance part of the potential may be affected by 
HYP. 



D. Results 

We obtain very good signals for the potentials and the 
forces defined by differentiating the potential with re- 
spect to R. The results at /3 = 2.35 on the 24 3 x 6 
lattice and on the 36 3 x 6 lattice, and at j3 = 2.43 on the 
24 3 x 8 lattice are plotted in Figs. [H [3] and 0] respectively. 
The q-q distances R of the potentials and the forces are 

improved to Ri and R = (^{G(R - a) - G(R)})~* , re- 
spectively. We fit these potentials to the function Vfn(R) 
in Eq. and extract the string tension and the Coulom- 
bic coefficient, which are summarized in Table Hill Since 
the potential and the force at = 2.20 on the 24 3 x 4 
lattice are already published in Ref. (3o| , we only present 
the fitting result for this data set. 
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FIG. 4: The same plot as in Fig.[2]at (5 = 2.43 on the 24 3 x 8 
lattice. 





Vna 


0.072(3) 


0.48(9) 


0.53(3) 


4.6 - 12.1 


1.03 




Va 


0.073(2) 


0.47(6) 


1.10(2) 


4.3 - 11.2 


1.03 




Vmon 


0.073(3) 


0.46(7) 


1.43(3) 


4.0 - 11.8 


1.01 




v ph 


-1.0(1) x 10~ 4 


0.0132(1) 


0.4770(2) 


6.4 - 11.5 


1.03 



TABLE III: Best fitted values of the string tension aa 2 , the 
Coulombic coefficient c, and the constant fia for the potentials 
Vna, Va, Vm 0n and Vph- 



24^ x 4 


aa 1 


c 


fj,a 


FR(R/a) 




Vna 


0.181(8) 


0.25(15) 


0.54(7) 


3.9 - 8.5 


1.00 


Va 


0.183(8) 


0.20(15) 


0.98(7) 


3.9 - 8.2 


1.00 


Vmon 


0.183(6) 


0.25(11) 


1.31(5) 


3.9 - 6.7 


0.98 


v ph 


-2(1) x 10" 4 


0.010(1) 


0.48(1) 


4.9 - 9.4 


1.02 


24 a x 6 


Vna 


0.072(3) 


0.49(6) 


0.53(3) 


4.0 - 9.0 


0.99 


Va 


0.073(4) 


0.41(7) 


1.09(3) 


3.7 - 10.9 


1.00 


Vmon 


0.073(4) 


0.44(10) 


1.41(4) 


3.9 - 9.3 


1.00 


v ph 


-1.7(3) x 10~ 4 


0.0131(1) 


0.4717(3) 


5.1 - 9.4 


0.99 



36 d x 6 



24 J x 8 



Vna 


0.0415(9) 


0.47(2) 


0.46(8) 


4.1 - 7.8 


0.99 


Va 


0.041(2) 


0.47(6) 


1.10(3) 


4.5 - 8.5 


1.00 


Vmon 


0.043(3) 


0.37(4) 


1.39(2) 


2.1 - 7.5 


0.99 


Vph 


-6.0(3) x 10~ 5 


0.0059(3) 


0.46649(6) 


7.7- 11.5 


1.02 



THE ABELIAN DUAL MEISSNER EFFECT 

Correlation function for the field profile around 
the q-q system 



Abelian dominance is seen again beautifully as in 
Sec. IIIII Moreover, we observe monopole dominance, 
i.e., the string tension of the static potential from the 
monopole Polyakov loop correlation function is identical 
to that of the non- Abelian static potential, while the po- 
tential from the photon Polyakov loop correlation func- 
tion contains no linear part. It is remarkable that Abelian 
dominance and monopole dominance for the string ten- 
sion are almost perfect as explicitly shown in Fig. [5j 
which also show the good scaling behavior with respect to 
the change of lattice spacing. We do not see the volume 
dependence of the string tension as shown in Fig. [6] 



These results suggest that although the lattice 
monopoles defined in Eq. ([5]) are gauge-dependent, they 
contain physical gauge-invariant pieces responsible for 
confinement, which show up after taking the gauge aver- 
age. 



We investigate the correlation function [43j, |44[ be- 
tween a Wilson loop W and a local Abelian operator O 
connected by a product of non-Abelian link variables 
(Schwinger line) L, 



(C(r)) 



(Tr [LW(i2,T)£t CT iO( r )]; 
(Tt[W(R,T)}) 



(13) 



A schematic figure is depicted in Fig. [7] 

We shall use the cylindrical coordinate (r, (f), z) to 
parametrize the q-q system, where the z axis corresponds 
to the q-q axis and r to the transverse distance as shown 
in Fig. [5] We are interested in the field profile as a func- 
tion of r on the mid-plane of the q-q system. 



B. Simulation parameters 

In this computation, we employ the improved Iwasaki 
gauge action [45| with the coupling constants (3 = 1.10 
and 1.28 on the 32 4 lattice, and f3 = 1.40 on the 40 4 
lattice in order to investigate the scaling behavior of the 
correlation functions with less finite lattice cutoff effects. 
Simulation parameters are listed in Table HVl The lattice 
spacings are determined so as to reproduce the physical 
string tension yfa = 440 [MeV] . To improve the signal- 
to-noise ratio, the APE smearing is applied to the spatial 
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FIG. 5: The o(/3) dependence of the square root of the non- 
Abelian, Abelian and monopole string tensions for the same 
temperature T = 0.8T C . The bottom axis for a set of three 
data points at the same lattice spacing is slightly shifted to 
distinguish each other. 
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FIG. 7: A schematic figure for the connected correlation func- 
tion. 
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FIG. 8: Definition of the cylindrical coordinate (r, <f), z) along 
the q-q axis. 



FIG. 6: The volume dependence of the square root of the 
string tensions on the 24 3 x 6 and 36 3 x 6 lattices at (3 = 2.35. 



TABLE IV: Simulation parameters for the measurement of 
the field profile, n and a are the number of smearing steps 
and the smearing parameter, which are optimized to obtain 
reasonable signals. 



p 


V 


a{0) [fm] 


Aconf 


n 


a 


1.10 


32 4 


0.1069(8) 


5000 


80 


0.2 


1.28 


32 4 


0.0635(5) 


6000 


80 


0.2 


1.40 


40 4 


0.0465(2) 


7996 


80 


0.2 



links of the Wilson loop [46j. We use the Wilson loop 
W(R = 3,T = 5) at f3 = 1.10, W(R = 5,T = 5) at 
(3 = 1.28 and W(R = 7,T = 7) at (3 = 1.40. Note that 
the physical q-q distance is the same (i? = 0.32 [fm]) for 
these Wilson loops. 



TABLE V: The parameter A corresponding to the penetration 
depth. 



p 


W(R,T) 


A [fm] 


Cl 


CO 


1.10 
1.28 
1.40 


W(3,5) 
W(5,5) 
W{7J) 


0.1075(13) 
0.1077(14) 
0.106(4) 


6.09(18) x 
1.024(14) x 10~ 2 
3.40(17) x 10 -3 


9(2) x 10 _b 
4.6(8) x 10" 6 
1.6(8) x 10" 5 



C. The penetration depth 

We measure all cylindrical components of the color- 
electric fields O(s) = EAi{s) — &4i(s). The results are 
plotted in Fig. [5] We find that only Ea z has correlation 
with the Wilson loop. We then fit {E\ z (r))w to a func- 
tion f(r) = ci exp(— r/A) + cq and find that the profile of 
{Ex z {f))w is well described by this functional form, i.e., 
the color-electric field is exponentially squeezed. The fit- 
ting curves are also plotted in Fig. [5J The parameter A 
corresponds to the penetration depth and the values for 
three gauge couplings are summarized in Table [V] and 
plotted in Fig. [TU] as a function of lattice spacing a((3). 
We find that the penetration depth A shows the good 
scaling behavior. 
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FIG. 9: The profile of the color-electric field E A at = 1.10 
(upper), /3 = 1.28 (middle) and P = 1.40 (lower). 



D. The dual Ampere law 

To see what squeezes the color-electric field, we study 
the Abelian (dual) Ampere law derived from the defini- 
tion of the monopole current in Eq. l[5]). 



Vx£ A = d 4 B A + 2nk 



(14) 



where Bm{s) = (l/2)eijkQjk(s)- The correlation of each 
term with the Wilson loop is evaluated on the same mid- 
plane of the q-q system as for the profile measurements 
of the color-electric field. We find that only the az- 
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FIG. 10: The penetration depth A as a function of lattice 
spacing a(P). 



imuthal components are non- vanishing, which are plotted 
in Fig.[TT] Note that if the color-electric field is purely of 
the Coulomb type, the curl of the electric field is zero. On 
the contrary, the curl of the electric field is non- vanishing 
and is reproduced mostly by the monopole currents. In 
any case, the dual Ampere law is satisfied, which is a clear 
signal of the Abelian dual Meissner effect. This result is 
quite the same as that observed in the MA gauge [13, [3 • 



E. The coherence length 

Let us estimate the coherence length by evaluating the 
correlation function between the squared monopole den- 
sity O(s) = k*(s) and the Wilson loop g2j. To measure 
such a correlation function, we use the disconnected cor- 
relation function, since the Schwinger lines are cancelled 
and the connected correlation functions for the squared 
monopole currents are automatically reduced to the dis- 
connected ones. Simulation parameters, the lattice vol- 
ume and the gauge couplings are the same as the mea- 
surements of the color-electric field profile, but the num- 
ber of gauge configurations is increased, N cou i = 5500 
for P = 1.10 and N con{ = 11887 for (3 = 1.40. For 
P = 1.28 we use the same number of configurations 
N con f = 6000. The physical q-q distance is again fixed to 
R = 0.32 [fm]. To reduce the noise, we further produce 
Nrqt = 100 gauge copies for each independent config- 
uration by applying the random gauge transformations 
and take gauge-averaging. 

The results are plotted in Fig. [T3J We then fit 
the profile of (k£(r))w to the functional form g(r) — 
c[ exp(— V2V/0 + Cg, where the parameter £ corresponds 
to the coherence length. We obtain the values for £ as 
summarized in Table IVT1 The coherence length shows the 
scaling behavior as demonstrated in Fig.[T3]as a function 
of lattice spacing a(P). 
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FIG. 11: Tests of the dual Ampere law at = 1.10 for W{R = 
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TABLE VI: The parameter f/v2 corresponding to the coher- 
ence length. 
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C/V2 [fm] 
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F. The vacuum type 



FIG. 12: The profile of the squared monopole currents at 
(3 = 1.10 (upper), 1.28 (middle) and 1.40 (lower). 



Taking the ratio of the penetration depth and the co- 
herence length, the GL parameter \/2/c = A/£ can be esti- 
mated, which characterizes the type of the superconduct- 
ing vacuum. The results are plotted in Fig. [TH against 
lattice spacing a 09). We obtain V2k = 1.04(7), 1.19(5) 
and 1.09(8) for (3 = 1.10, 1.28 and 1.40, respectively. 

We find that the GL parameter shows the scaling be- 
havior and the value is about one. This means that the 
vacuum type is near the border between the type 1 and 2 
dual superconductor. However, we note that the physical 
spatial size of the Wilson loop used in the present sim- 
ulations is still small (i? = 0.32 [fm]). Clearly, further 
quantitative studies with larger Wilson loops are needed 



to determine the definite value. 



VI. NON-ABELIAN COLOR CONFINEMENT 

Let us consider what is induced from the above numer- 
ical results. 

Since gauge fixing is not applied in these computa- 
tions, Abelian fields in any color directions are equivalent. 
Thus, our result is interpreted as that the color-electric 
fields in all color directions are squeezed and the Abelian 
(monopole) string tensions in all color directions are the 
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FIG. 13: The coherence length £ as a function of the lattice 
spacing a(/3). 
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FIG. 14: The GL parameters as a function of the lattice spac- 
ing a(J3). 



same as the non-Abelian string tension. This indicates 
that QCD contains a gauge- invariant Abelian mechanism 
of confinement which is not related to the specific gauge 
fixing. Namely Abelian monopolcs in three color direc- 
tions arc condensed in the vacuum of the confinement 
phase of 577(2) QCD. 

Let us denote quark fields having charge 1/2 and —1/2 
in the a 3 direction, respectively, as u 3 and d 3 . Then local 
mesonic states, u 3 u 3 and ^3^3, are Abelian color neutral 
in the 173 direction. Consider next 



u 3 + d 3 u 3 - d 3 

iti = ■=— , di = 



u 2 



V2 
iu 3 + d 3 

~7T 



d 2 = 



V2 ' 
iu 3 - d 3 

V2 • 



U\ and d\ (1*2 and d 2 ) are quark fields having charge 
1/2 and —1/2 in the o\ (02/2) direction. Using these 



expressions, the quark-gluon coupling term is written as 
= z\{u 3lll d 3 + d 3 ^u 3 )Al 
~i\{u 3 l^d 3 - d 3 jf,u 3 )Al 
+ r\(u 3 j„u 3 - d 37fl d 3 )Al (15) 



2^17^1 



+^(u2lf,u 2 - d27f J ,d 2 )Al 

+ ^0"37p"3 - d 3 j^d 3 )Al, (16) 



where the first equation (j 1 5[) is expressed in terms of u 3 
and d 3 alone. Consider local mesonic states U\U\ and 
d\di {112U2 and d 2 d 2 ) which are Abelian color neutral in 
the u 1 ((72 ) direction. When we look at the states %\U\ 
and d\d\ in the a 3 direction, they are written as the sum 
of color- neutral and color-charged states: 



U\U\ 



■^(u 3 u 3 + d 3 d 3 + u 3 d 3 + d 3 u 3 ) , (17) 



d idi = -^{u 3 u 3 + d 3 d 3 ~ u 3 d 3 - d 3 u 3 ) . (18) 

The same observation applies to the color-neutral states 
U2U2 and g?2<^2 in the 02 direction. However, we find that 



uiui + d\d\ = U2U2 + d 2 d 2 = u 3 u 3 + d 3 d 3 



(19) 



are Abelian color neutral in all color directions. The state 
(Un]) is nothing but the non-Abelian color singlet state. 

This example tells us that the Abelian color-neutral 
state in any color directions corresponds to the physi- 
cal non-Abelian color-singlet state. Hence, the confine- 
ment of non-Abelian color charges can be explained in 
terms of the Abelian dual Meissner effect due to Abelian 
monopoles. To the authors knowledge, this is the first 
paper that explains the confinement of non-Abelian color 
charges only in terms of the Abelian dual Meissner effect. 



VII. CONCLUDING REMARKS 

We make some concluding remarks. The Abelian 
gauge fields extracted from the thcrmalized non-Abelian 
link fields contain originally topological monopoles re- 
sponsible for the confinement mechanism of non-Abelian 
color charges even in the continuum limit. Our results 
presented in this paper are almost the same as those ob- 
tained in the maximally Abelian gauge. This suggests 
that the MA gauge fixing is the easiest method to ex- 
tract the physical ingredients of the monopoles, since we 
do not need very precise time-consuming simulations in 
the MA gauge as done here. 

In the lattice Landau gauge, it is known that no 
monopoles exist [48| if one uses DeGrand-Toussaint def- 
inition and the magnetic displacement current takes a 
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role of monopoles in the dual Ampere law. How to inter- 
pret the existence of a gauge-invariant Abelian confine- 
ment mechanism in the framework of the Landau gauge ? 
Abelian monopoles are as a whole gauge-variant with- 
out gauge fixing, but they may contain a gauge-invariant 
physical component and a gauge-variant unphysical one. 
The compatible interpretation would be that the Lan- 
dau gauge is a special gauge in which the unphysical 
gauge variant piece apparently cancels the physical one 
in the DeGrand-Toussaint monopoles, but at the same 
time, the role of physical monopoles are carried by the 
color-magnetic displacement current, which is just a mat- 
ter of definition of monopoles on the lattice. On the other 
hand, in the MA gauge, the main part of the DeGrand- 
Toussaint monopoles is a physical component. 

If there exist physical gauge-invariant ingredients of 
Abelian monopoles, one could observe them in the real 
experiment [491 ]. To find the effect in the confinement 



and also in the deconfinement phases is a very interesting 
topic in the future. 
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